Measurement and control of spatial qubits generated by passing 

photons through double-slits 

Gen TaguchiJ*| Tatsuo Dougakiuchi, Nobuaki Yoshimoto, Katsuya 
Kasai, Masataka Iinuma, Holger F. Hofmann, and Yutaka Kadoya 
Graduate School of Advanced Sciences of Matter, 
Hiroshima University, Higashihiroshima, 739-8530, Japan 
(Dated: May 8, 2008) 

Abstract 

We present an experimental study of the non-classical correlations of a pair of spatial qubits 
formed by passing two down-converted photons through a pair of double slits. After confirming 
the entanglement generated in our setup by quantum tomography using separate measurements 
of the slit images and the interference patterns, we show that the complete Hilbert space of the 
spatial qubits can be accessed by measurements performed in a single plane between the image 
plane and the focal plane of a lens. Specifically, it is possible to obtain both the which-path and 
the interference information needed for quantum tomography in a single scan of the transversal 
distribution of photon coincidences. Since this method can easily be extended to multi-dimensional 
systems, it may be a valuable tool in the application of spatial qudits to quantum information 
processes. 
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I. INTRODUCTION 



Quantum information science makes use of the entanglement between well-defined two 
dimensional qubits or d-dimensional qudits to perform tasks that could not be performed 
by a corresponding classical system. In optical implementations of quantum information 
technologies, a commonly used source of entanglement is the generation of photon pairs by 
spontaneous parametric down-conversion (SPDC) [l, 2|. 

In many applications, photon polarization is used to define a qubit, since it provides a 
natural two level system that is comparatively easy to control by using birefringent optical 
elements. However, photon pairsproduced by SPDC are also entangled in their spatial 
and temporal degrees of freedom . Since these degrees of freedom are naturally 

continuous, it is necessary to introduce additional constraints in order to define a qubit or 
qudit system. The advantage of this approach is that it is fairly easy to extend the method 
;o higher dimensions, a possibility that may simplify some quantum information processes 

3, y,y. 

One widely used method of defining spatial qudits by discretizing the transversal degrees 
of freedom is the selection of angular momentum eigenstates corresponding to photons in 
Gauss-Laguerre modes 1(3, 11, 12, Q, Q- However, it is in principle not necessary to base 



the selection of modes on symmetries, since the spatial entanglement applies to arbitrary 
selections of orthogonal modes. It is therefore possible to define qudits by simply selecting 
a sufficiently narrow spatial "pixel" for each basis state of the qudit [l^|. Alternatively, it is 
possible to concentrate on only one spatial dimension. As was demonstrated by Neves and 
co-workers, entangled spatial qudits can then be obtained by using the familiar slit arrays 
used in basic demonstrations of optical interference 161 ] . 

An interesting and important aspect of the use of multi-slits to define spatial qudits is the 
fact that the photons continue to propagate in continuous free space after their wavefunction 
has been projected into a two dimensional Hilbert space [l_7j] - This means that the spatial 
qudits continue to exist in the infinite dimensional Hilbert space of transverse position and 
momentum. Specifically, the only directly observable physical property of the qudits is their 
transverse position, so that the measurement and control of the qudit states has to be based 
on the evolution of their spatial wavefunction. It is therefore of great interest to explore 
the possibilities of preparing and characterizing different multi-slit qudit states using their 
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propagation in space. 

In this paper, we present a thorough investigation of the entanglement between a pair of 
double-slit qubits based on measurements of their correlated spatial patterns. In particular, 
we show that measurements performed between the focal and image planes of a lens simulta- 
neously provide which-path information that distinguishes between the slits and interference 
information that identifies the quantum coherence between the slits. It is therefore possible 
to scan the whole surface of the qubit Bloch sphere by varying the detector position in 
this intermediate plane. We can use this method to prepare arbitrary superposition states, 
and to perform quantum tomography of the states thus prepared by scanning only a single 
transverse pattern. 

The rest of the paper is organized as follows. In section [III we present the experimental 
setup and report the results of quantum tomography performed by measuring count rates 
corresponding to eigenstates of the Pauli operators in the image and focal planes. In section 
IIII[ we analyze the effects of a measurement between the focal and image planes on the spatial 
qubit and show how it can be used for qubit preparation and for single-scan tomography. 
In section HVl we present the experimental results of single-scan tomography for conditional 
states prepared by measurements in the other arm and reconstruct the density matrix of the 
entangled state once more from these results. The results are compared with those from the 
Pauli operator measurements, and possible experimental problems are discussed. Section [V] 
concludes the paper. 



II. GENERATION OF ENTANGLED SPATIAL QUBITS 

In this work, we use double-slits to define the spatial qubits. It is convenient to express 
the qubit states in the {\l), \r)} basis corresponding to photons passing through the left or 
right slits as shown in fig. [TJ In this basis, the entangled state ideally generated by our 
setup is given by 

|^ s iit s ) = ^(|0A|r) B + |r) A |0 B ), (1) 

where the suffixes A and B denote two photons found in different arms of our setup. In 
this section, we explain the setup in detail and present an experimental confirmation of the 
entanglement of our source by quantum tomography. 
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FIG. 1: Definition of the entangled qubits. The slits in each arm A, B define the basis states 
{\l), \r)} corresponding to the slit a photon passes through. 
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FIG. 2: Experimental setup. Note that basis states in the B arm are exchanged due to the reflection 
at the PBS. The detector system is composed of the parts shown in the inset. 

A. Experimental setup 

Figure [2] shows our experimental setup. The entangled photons were generated in type 
II SPDC by a 405 nm pump beam from a 45 mW CW laser incident on a 5 mm-thick /3- 
Barium Borate (BBO) crystal. The BBO crystal was placed in the collinear condition. The 
photons were separated into two different directions by a polarizing beam splitter (PBS). In 
each arm, a double-slit was placed in the focal plane of a lens in order to select a pair of 
transverse momenta of the photons. The slit width was 40/xm and the distance between the 
slits was 150^tm . A lens of focal length / = 50 mm was placed at a distance of 2f from the 
double-slit. A detector system was placed in a plane at a distance of z from the lens. 
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FIG. 3: Measurements of Pauli operators, (a) a z can be measured in the image plane, (b) Any 
linear combination of a x and a y can be measured in the focal plane, (c) Detectors in the image 
plane and the focal plane measure the states on the poles and the equator of the Bloch sphere. 

The detector system consists of a single slit with a width of 40/im, a band pass filter (810 
nm, band width of 10 nm), an objective lens to couple the photons into a multimode fiber, 
and a photon detector (Perkin-Elmer SPCM-AQR-14) as shown in the inset of fig. [2j The 
coincidence counts between the detectors in the two arms were recorded. 

B. Measurement of Pauli operator averages in the focal and image planes 

The spatial qubits can be characterized in terms of the Pauli operators, a x , a y and 
o z . In the {|/), |r)} basis, they are a x = \l)(r\ + \r)(l\, cr y = — i(\l)(r\ — \r)(l\), and o z = 
— \r)(r\. To measure a z , we have to distinguish the two slits. Since the image of the 
double-slit appears in the image plane at a distance of z = 2f from the lens (fig. [3^a)), we 
can do this by placing the detectors at the positions corresponding to each slit in the image 
plane. To measure a x and a y , we have to obtain the phase information of the interference 
pattern of the double slits. Since the interference patterns are observed in the focal plane 
at a distance of z = f from the lens (fig. E](b)), we can do this by placing the detectors 
at appropriate positions in the focal plane. E. g. for a photon in the state \l) + |r), an 
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interference pattern appears as indicated in fig. [3](b). The a x operator averages can then 
be estimated from the difference between the count rates at the center and at the first node 
of this interference pattern. These two detector positions thus correspond to projective 
measurements of the o x eigenstates, \l) + |r) and |Z) — |r). Likewise, a measurement of 
the a y eigenstates can be realized by placing the detectors at the positions corresponding to 
\l)+i\r) and \l)—i\r). We can therefore realize the Pauli operator measurements by detecting 
photons at each of the six positions in the two planes corresponding to the eigenstates of 
a x , <T y , and a z , as shown in fig. [3](c). 

C. Results of quantum tomography 

To measure the correlations between the spatial qubits, the detector in arm B was fixed 
at a position corresponding to one of the eigenstates of the Pauli operators. In arm A, the 
detector position x was scanned in either the focal plane or the image plane. First, the 
scanning detector was placed in the focal plane. When the fixed detector was placed at the 
appropriate position for a measurement of a a x or a y eigenstate, a conditional interference 
pattern was observed. As the fixed detector position was moved, the interference pattern 
changed its phase without losing visibility. The interference patterns for the four eigenstates 
of a x or o y are shown in fig. IH(a) and (b). When the fixed detector was set at a position 
corresponding to a a z eigenstate, the interference pattern disappeared as shown in fig. 11(c), 
that is, we see no correlation between <j z and a x ja y . Next, the scanning detector was placed 
in the image plane. The results are shown in fig. [5j The o z -o z correlation was clearly 
observed, whereas no u z -u x and u z -u y correlation was found, as expected. 

The density matrix can be reconstructed by using the data extracted from the correlation 
experiments and listed in Table [1(a). Specifically, the data for the a x measurement in the 
scanning arm A was obtained from the scanning detector positions at or 135/im from the 
center in fig. HI and the data for the o y measurement in the scanning arm A was obtained 
from the scanning detector positions at ±67/im from the center in fig. HI The data for the 
a z measurement was taken from the centers of the slit images at ±70/im in fig. [51 To obtain 
the probabilities, it is necessary to normalize the results. Since the coincidence count rates 
decrease as the detectors move away from the optical axis, it is also necessary to compensate 
the resulting detection efficiency differences between the \l) + |r) and \l) — \r) states using the 
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Detector Position (|im) 



FIG. 4: Coincidence count rates as a function of scanning detector position. The dots and squares 
show the experimental data and the solid curves are guides for the eyes. The detector positions 
at and 135//m correspond to a a x measurement in the scanning arm A. The positions ±67/^m 
correspond to a a y measurement in the scanning arm A. (a) Interference correlation. The detector 
in arm B was set to the a x eigenstates \l) + |r) (dots) and |/) — |r) (squares), (b) The detector 
in arm B was set to the a y eigenstates \l) + i\r) (dots) and \l) — i\r) (squares), (c) Interference - 
which-path correlation. The detector in arm B was set to a position corresponding to one of the 
o z eigenstates 



theoretically expected ratios between the corresponding peaks in fig. H](a). The normalized 
results are shown in table [Jb). From these probabilities, we can reconstruct the density 
matrix. In the {\ll), \lr), \rl), \rr)} basis, it is given by 

/ 0.003 -0.005 - 0.007z -0.006 + 0.000* 0.002 - 0.006z \ 

-0.005 + 0.007z 0.498 0.463 - 0.024i 0.009 + O.OOlz 

-0.0060.000i 0.463 + 0.024i 0.497 0.008 - 0.007i 

\ 0.002 + 0.006i 0.009 - O.OOli 0.008 + 0.007z 0.002 / 



Pab 



(2) 
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FIG. 5: Coincidence count rates as a function of scanning detector position. The dots and squares 
show the experimental data and the solid curves are guides for the eyes, (a) Which-path correlation 
measurement. The detector in arm B was set to the a z eigenstates \l) (dots) and \r) (squares), (b) 
One of the which-path - interference correlation measurements. The detector in arm B was set to 
a position corresponding to one of the a x eigenstates + |r)). 




FIG. 6: Real and imaginary part of the reconstructed density matrix. 

Figure [6] shows an illustration of the density matrix elements. The fidelity of the ideal 
state given in eq. (pQ) is F = {\J/ S ii ts |pa_b | ^siits) =0.961, indicating that we have successfully 
generated an entangled state close to the theoretically expected one. 
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TABLE I: Coincidence counts. The columns give the positions of the scanning arm A and the lines 
give the setting of the fixed arm B settings. 

(a) Raw data. 
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595 


290 


313 
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\l)-i\r) 
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(b) Normalized data 





10 


r) 


10 + 1*-) 


10 - |r) 
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\l)-i\r) 


10 


0.003 


0.497 


0.253 


0.262 


0.252 


0.248 


r) 


0.498 


0.002 


0.252 


0.233 


0.251 


0.249 


\l) + \r) 


0.245 


0.255 


0.486 


0.017 


0.275 


0.227 


10 - |r> 


0.258 


0.242 


0.017 


0.480 


0.243 


0.255 


\l)+i\r) 


0.258 


0.256 


0.254 


0.262 


0.484 


0.025 


\l)-i\r) 


0.238 


0.248 


0.256 


0.228 


0.014 


0.477 



III. MEASUREMENT BETWEEN THE FOCAL AND IMAGE PLANES 



The measurement of an arbitrary superposition state of the spatial qubit may play an 
important role in the application of the spatial qubit entanglement. In this section, we show 
that such measurements can be realized by placing the detector between the focal and image 
planes, as shown in fig. [7J The positive operator valued measure describing the measurement 
in this intermediate plane is derived and its application to state preparation and quantum 
tomography is discussed. 
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FIG. 7: Simultaneous measurement of interference and which-path information realized by placing 
the detector between the focal and image planes. 

A. Positive operator valued measure in the double-slit basis 

The quantum states of the photons at the double-slits can be expressed in terms of a 
transverse spatial wavefunction (p(0;x). If we place a lens of focal length / at a distance L 
from the slits, then the transverse wavefunction ^p(L; x) in the plane in front of the lens can 
be calculated using the Fresnel-Kirchhoff diffraction integral 

<p(L;x) = JJ-e i ™e i ™ a? [ ^x'Y^e-^'dx'. (3) 
V AL J So 

where the integral is performed over the one-dimensional slit plane. The diffraction at 
the lens modifies the free space propagation so that the wavefunction at the distance z 
from the lens corresponds to the wavefunction at an effective propagation length of R = 
(Lf + zf — Lz)/(z — f) from the double-slits, reduced in size by a factor of (z — /)//. Up 
to an arbitrary phase factor, the wavefunction at z is then given by 



ip(L + z;x) = \l XR ^_ f) e~ l ^ x2 [ <p(0; x^e'^' 2 e'^^' dx'. (4) 



So 



The wavefunction of a slit of width a at a distance of r n from the optical axis is given by 



y n {L + z; x) = A / -—L—e- 1 ^* 2 f ^ VTO* ( *' +r "W, (5) 

V ^Ra{Z - /) J -a/2 

where we have assumed that the wavefunction in the slit is uniform. If the detector plane 
is far enough from the image plane, so that R > a 2 /X, eq. (jHJ) can be simplified to the 
conventional slit diffraction pattern given by the sine function sinc(x) = sinx/x, 



[K ( 2r n \ 
(p n (L + z;x) = y — exp I — i Kx J sine 



K ( x + Z - - ' r r 



(6) 
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where K = iraf/XR(z — /) defines the scale of the single slit diffraction pattern. 

For the double-slit system, r\ = —d/2 and r> = d/2. Each basis state of the qubit system 
thus corresponds to a well-defined wavefunction in the plane at z from the lens. These two 
wavefunctions define the two dimensional Hilbert space of the qubit in the measurement 
plane. It is then possible to express the effect of a measurement of x in this plane by a 
projection onto a non-normalized state in the Hilbert space of the qubit, 

\m{x)) = tf{L + z;x)\l)+ip* r {L + z;x)\r). (7) 

We can therefore identify each measurement result x with a specific projection state of the 
spatial qubit, corresponding to a point on the Bloch sphere. The positive operator valued 
measure describing the qubit measurement is given by the measurement operators 

M(x) = \m(x))(m{x)\. (8) 

For a qubit state given by a density operator p, the probability of photon detection at the 
position x is 

P(x)=Tr[Mp] 7 (9) 

where the normalization of the wavefunctions (pi(L + z\ x) and (p r (L + z\ x) ensures that the 
set of operators {M(x)} satisfies the completeness relation 

/oo 
M(x)dx = 1. (10) 
-oo 

We can easily extend this method to the analysis of mult i- dimensional qudits. All we 
have to do is increase of the number of slits to N. The projection state of eq. (J7|) is then 
given by a superposition of all slit states, 

N 

\m^(x)) = Y^<Pn(L + z ] x)\n) t (11) 

n=l 

where ip n (L + z; x) represents the wavefunction of a photon originating from the nth slit, as 
given in eq. (jSJ). 

B. Accessibility of arbitrary quantum states 

In our setup, L is equal to 2f. In this case, the effective propagation length is R = 
f(2f — z)/(z — /), so we can access any value of R between R = (image plane) and 
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(a) 3 (b) 2 (c) 




z=\.lf z=\.6f z=\.9f 

R=9f R = 2f/3 R = f/9 

FIG. 8: Measurement states \m{x)) on the Bloch sphere. The calculation was performed with the 
distance z of the detector from the lens fixed, while the transverse position x was varied until the 
intensity from one of the slits went down to almost zero, (a)-(c) show the trajectories of the Bloch 
vectors for the respective distance z of the detector from the lens. 

R —>■ oo (focal plane) by placing the detectors at a distance z from the lens between z — f 
and z=2/. For each detector setting, the projection states \m(x)) given by eq. ([7]) trace out 
a trajectory on the Bloch sphere as x is scanned. The trajectories for three specific detector 
positions z are shown in fig. [5J In each case, the Bloch vector moves from the upper 
hemisphere to the lower hemisphere rotating around the a z axis as the detector position x is 
scanned. For z close to /(fig. Ufa)), the Bloch vectors are concentrated around the equator. 
They spread out as z increases and accumulate at the poles as z approaches 2/ (fig. [HJc)). 
In the intermediate case shown in fig. [SJb), the Bloch vectors are almost equally distributed 
over the Bloch sphere, indicating that any state on the Bloch sphere can be measured by an 
appropriate choice of detector position (z, x). 

Using our source of entangled qubits, we can prepare an arbitrary state of one spatial 
qubit by post-selecting a specific measurement outcome of the other qubit. For the ideal 
state of eq. ([I]), the non- normalized state prepared in A by post-selecting a measurement 
result of x in B is 

|*prcp)A = B{m{x B )\%> ' s i it } AB = (if r (L + Z)X B )\l) + <fl(L + Z\X B )\t)) ■ (12) 

The Bloch vector of this state is the mirror image of the Bloch vector of \m(x)) at the 
X — Y plane. The trajectory of states that can be prepared in A by measurements in B are 
therefore equivalent to the trajectories shown in fig. [HI 
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C. Spatial patterns of density matrix elements 



From the count rate pattern obtained by scanning the detector position x in a plane 
between the focal and image planes, we can know the state of the spatial qubit. In the 
double-slit basis, the density matrix of the spatial qubit at the double-slit is given by 




P 



For a detector position z, the positive operator valued measure of eq. 
of detection at a point 



(13) 

gives the probability 



P(x) 



Tr 



pM 



pu \ipi(L + z;x)\ 2 + pi r <p*(L + z;x)ip r (L + z;x) 
+p r i tpi(L + z; x)tp*(L + z; x) + p rr \<f r (L + z; x)\ 2 . 



(14) 



Each density matrix element is thus connected to a unique pattern defined by (pi and (p r 
Using eq. we obtain 

2 



P(x) 



K 



Pa- 



li 



sine ( Kx — 



2d 



K 




+p ir —exp 


7T 




K 




+p r i— exp 


7T 




K 




~\~ Prr 

7T 


sine 



—i — Kx ) sine ( Kx 

a 

' 2d 

i — Kx } sine I Kx 



A 4) 



A(f) 



sine ( Kx 



Ao\ . ... 

sine \ K x + 

2 J V 2 



Ad) 
A0 



(15) 



where A<p = (z — f)Kd/ f describes the spatial displacement of the diffraction patterns of the 
two slits, that is, A<p is proportional to the ratio of the separation between the slit images 
and the width of the diffraction pattern of a slit. The closer the detection plane comes to 
the focal plane, the smaller the separation between the slit images becomes. 

In general, the detection probability P(x) is a linear function of the four density matrix 
elements. Therefore, it is in principle always possible to invert and solve the equation for the 
elements, so that the density matrix elements can be determined from specific integrals of 
P(x). In the case of the qubit in eq. ( fTBT) . the complete density matrix can be determined by 
only two coefficients, the difference between the diagonal elements pu — p rr and the complex 
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sine 2 ( Kx — | — sine 2 ( Kx + 



P(x) dx 



off-diagonal element pi r . In P(x), pu — p rr appears as a difference between the squared sine 
functions describing the diffraction patterns of the two slits and p\ r appears as an oscillating 
interference pattern with an envelope given by the product of the two sine functions. If the 
oscillation is sufficiently fast {d ^> a), the inversion can be performed by integrating the 
product of P(x) and the patterns corresponding to pu — p rr and pi r , respectively. The result 
of the inversion is then given by 

3 

PU-Prr ~ 2(1 _ /3) J L 2 J ^— ■ 2 

3 f f 2d \ ( A<A ( A<A . , 
Pi r = — J exp I i — Kx J sine I Kx — \ sine I Kx H — — \ P[x) dx, 

3 

where f3 = — — (1 - sinc(2A</>)). (16) 
2A(p 2 

Since pu + p„ = 1, these two patterns completely define the density matrix of the qubit. For 
higher dimensional qudit systems, similar relations can be derived for all iV 2 density matrix 
elements of an iV-slit system. Thus the density matrix of the spatial qudits can be directly 
identified with the corresponding spatial patterns observed in the measurement. 



IV. EXPERIMENTAL RESULTS 



In this section, we present the experimental results of the measurements with the detector 
placed between the focal and image planes. Conditional measurements of the pattern P(x) 
in arm A were performed with a fixed detector position x in arm B. From this data, we 
can determine the conditional density matrices of the states in A. The density matrix of the 
complete two qubit system is then reconstructed from all 6 conditional measurements and 
the result is compared with that obtained in section [Til 



A. State preparation and conditional density matrices 

The experimental setup was as described in section [TTJ except that the slit width in the 
detector system was 20pm. The distance z was set to 1.8/ in each arm. Six different detector 
positions were chosen for arm B such that the corresponding Bloch vectors are close to a 
regular octahedron, as shown in fig. [91 The detector in A was then scanned to obtain the 
coincidence count data shown in fig. [TUl The distributions measured in each scan correspond 
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FIG. 9: Six Bloch vectors corresponding to the detector positions in the fixed arm B. 
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FIG. 10: Coincidence counts data. Each measurement was obtained by placing the fixed detector 
in a different position. The insets show the corresponding conditional density matrices. 



to a one qubit density matrix, where the diagonal elements determine the bias between the 
right and left side of the distributions and the off-diagonal elements determine the visibility 
and the phase of the interference patterns. 

We can understand the results of the conditional measurements from the viewpoint of 
quantum state preparation as explained in section IIIIB1 In our experiments, the detector 
position in the fixed arm B selects the prepared state in the scanning arm A. We can 
check how well the conditional state in arm A corresponds to the expected one described 
by eq. ( fl2|) by reconstructing the conditional density matrix and determining the fidelity. 
To reconstruct the density matrix and to determine the actual value of the experimental 
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TABLE II: Fidelities of expected prepared states and experimentally obtained states in fig. [TU] 



case 


(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


fidelity 


0.887 


0.861 


0.841 


0.841 


0.871 


0.912 



parameter L at the same time, we performed a least square fit based on eq. f|T4|) . The result 
of the fits and the corresponding density matrices are shown in fig. [TUJ The fidelities of the 
expected prepared states given by eq. ( fl~2l) are given in table [III The results show that the 
states prepared in A by measurements in B are close to the intended ones. 



B. Reconstruction of the density matrix from conditional scans 



The non-normalized conditional density matrices pa{x%) obtained by setting the detector 
in arm B to Xi depend on the density matrix pab of the entangled qubit according to 

PA{xi) = B(m(xi)\pAB\m(xi)) B - (17) 

In principle, four different measurement points Xi are enough to reconstruct the density 
matrix in B. The entangled density matrix of A and B can be reconstructed from the 
conditional density matrices pa{xi) with the same linear coefficients used to reconstruct a 
density matrix in B from probabilities P(xi). Specifically, the reconstruction of ps from 
P{xi) and the reconstruction of pab from p^(:rj) are related by a set of reconstruction 
operators Aj with 

p B = Y^ p (?i% ( 18 ) 



Pab 



A,:. 



(19) 



Since we have chosen 6 measurement points in B, we have sufficient information to recon- 
struct the complete two qubit density matrix pab from the conditional density matrices 
Pa{xi). The result in the {|//), \lr), \rl), \rr)} basis reads 

/ 0.008 0.008 -0.012i 0.015 + 0.021a -0.018 - O.OOli \ 
0.008 + 0.012z 0.485 0.347 - 0.038i 0.002 - 0.027i 
0.015-0.021i 0.347 + 0.038i 0.469 0.008 + 0.005i 

\ -0.018 + O.OOli 0.002 + 0.027i 0.008 - 0.0050z 0.038 
The fidelity of the ideal state given in eq. ([T|) is F — (^sUtslPAsl^siits) = 0.824. 



Pab 



(20) 



/ 
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C. Discussion of the results 

At present, the agreement with the theoretical prediction 1 is not as good as that of 
the Pauli operator measurements described in section [III It should therefore be possible 
to improve the experiment by identifying the sources of the additional errors. Since there 
are some discrepancies between the experimental data and the theoretical fit in fig. [10l one 
problem may be that the experimental parameters used in the theory need to be adjusted. 
Moreover, the visibilities may be underestimated by the fit as indicated by data points above 
and below the maxima and minima of the interference patterns. An increase in the spatial 
resolution of the measurement might solve this problem. 

The main merit of our method is that we can see both the which-path and the interference 
information in a single scan. Our experimental results clearly confirm that the coincidence 
count distribution contains all the data necessary to reconstruct the density matrix. This 
kind of single scan tomography thus provides a direct characterization of spatial qudits in 
terms of their actual physical properties in space. 

V. CONCLUSIONS 

We have presented a thorough experimental investigation of the entanglement between 
a pair of spatial qubits generated by passing down-converted photon pairs through double- 
slits. Pauli operator measurements indicate that we have achieved a fidelity of F — 0.961 for 
the maximally entangled two qubit states. We have then explored the possibility of accessing 
arbitrary superposition states by measurements between the focal and image planes. We 
can use the entanglement to prepare arbitrary states in arm A by measurements in arm B. 
The results have been verified by single scan tomography, where we identified each density 
matrix element with its distinct pattern in the measurement distribution. Fidelities between 
0.841 and 0.912 have been obtained for the conditionally prepared states. By combining the 
conditional results, we obtained a two qubit density matrix with a fidelity of 0.824. The 
decrease in the fidelity compared to the Pauli operator measurements suggests that further 
experimental improvements may be possible. 

The method introduced here allows us to access the full Hilbert space of the spatial qubits. 
It may therefore have applications in the realization of quantum information protocols using 
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double-slit qubits. Moreover, the extension of our method to iV-slit qudits is straightforward, 
opening the way to quantum operations in higher dimensional Hilbert spaces. 
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